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Abstract 

The decay modes B — > D^tP t are sensitive to charged scalar effects, such as 
the charged Higgs effects. In this paper we suggest a method to determine their 
effects by using the ratio of branching fractions and forward-backward asymme- 
tries. In particular, forward-backward asymmetries on B — > D^t(— > 7r^ r )z/ T , 
B D^t(-> pv T )v T , and B L>(*V(^ 

aiv T )v T play an important role, which 
discriminate the Standard Model from other New Physics scenarios. 
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1 Introduction 



Despite the fact that the Standard Model (SM) has been very successful in describing 
most of elementary particles phenomenology, the Higgs sector of the theory remains 
unknown so far, and there is no fundamental reason to assume that the Higgs sector 
must be minimal, i.e., only one Higgs doublet. The simplest extension compatible with 
the gauge invariance is called the two Higgs doublet model (2HDM), which consists 
of adding a second Higgs doublet with the same quantum numbers as the first one. 
Similarly, the minimal supersymmetric Standard model (MSSM) consists of adding a 
second Higgs doublet. In the MSSM, two Higgs doublets are introduced in order to 
cancel the anomaly and to give the fermions masses. The introduction of a second 
Higgs doublet inevitably means that a charged Higgs boson is in the physical spectra. 
So, it is very important to study effects of charged scalar particles. 

The branching fractions of B — > D£u e and B — > D*£vg have been measured in B 
Factories, where £ denotes e, /i or r. We define R(D^) as the ratios of the branching 
fractions, that is, 



R(D 



w = BjB -» D^TUr) 

' B(B -> D«(e or /i)v) 



Using the ratio of two branching fractions lowers the hadronic uncertainty. The theoret- 
ical predictions in the Standard Model using the heavy-quark effective theory(HQET) 
on5-> D^t&t are evaluated as 

#(£>)hqet = 0.310 ± 0.011, (2) 
R(D*) HQET = 0.253 ± 0.003. (3) 

These are consistent with the results in Refs. [HI2]- R{D) is also evaluated by using 
hadronic form factors computed in unquenched lattice QCD as R(D)\ at = 0.316(12)(7), 
where the errors are statistical and total systematic, respectively [3j. In Ref. [I], R{D) is 
evaluated by using results of HQET and lattice QCD as -R(-D)HQET+iat = 0.31(2). These 
theoretical predictions are consistent with each other within their errors. The recent 
experimental results of R(D^) by BABAR [5] are 

R(D) cxp = 0.440 ± 0.058 ± 0.042, (4) 
R(D*) cxp = 0.332 ± 0.024 ± 0.018, (5) 

which exceed the Standard Model expectations by 1.9a and 2.6a, respectively. 
In this paper, we consider an effective weak Hamiltonian such as 

% ^ t ) = A^[Oy L + m e C SR Os R + m t Cs L SL ] + H.c, (6) 

Vl = {afPMh^PLVtl (7) 
Os R = (cP R b)(£P L u e ), (8) 
Sl = (cP L b)(£P L u e ), (9) 

where are projection operators on states of positive and negative chirality. We 
assume that the neutrino helicity is only negative. This type or a more general one 



1 




Fi gure 1: @dW M 1S the angle between the direction of and M(~ n, p, ai) in the r — v T rest 
frame. 



has been studied in Refs. [2VHII6HT3] by using some observables, e.g, R(D^) and q 2 
distributions of R ratios and angular asymmetry on B — > D^*'tT> t , where q 2 = (j>b — 
Pd(*)?- 

Since a tauon decays into a light meson(lepton) with nutrino(s), the measurements 
of angular distribution for tauon on B — > D^*>rv T are difficult. However, angular de- 
pendence 011B4 D(*)tu t is important to search for the NP effect. So, we study rela- 
tions among the coefficients Cs RL and forward-backward asymmetries on B — > D^t(—> 
Txu r )u T , B — > D(*>t(-^ pv T )i> T , and B — >■ D^*'r(-^ aiv T )P T , and show that it is possible 
to determine them almost completely by using the ratios of the branching fractions and 
forward-backward asymmetries on these modes. 

2 R(D^) and Forward- Backward Asymmetries 

We use the quantities the ratios R(D^) defined as ([T]) and the forward-backward asym- 
metries A FB defined as 



(fo-f°i) dcos9 D( , ) M 



dr D(*),M 



A FB (D ( *\M) = ± ^L, (10) 

J- ; 



)(*) M\ - _ L . D( } - M 

d(*Km 

dT DW>M = dT{B £>«r(^ Mz/ r )z7 r ), (11) 
M = n, p or ai, (12) 



where 9 D (*) M is the angle between the direction of the and the M in the r — v T rest 
frame, as seen in Figured! The q 2 distribution and angular distribution on B — > D^rV T 
have been analyzed [21IHIHQ] - We can check the differential decay rate on B — > Z)r(-> 
7Tu T )u T in Ref [7]. The differential decay rates are written as 

dT(B-+DV>TU T ) = -!—d* 3 x l^.)^' 008 ^)! 2 ' (13) 

mB A T (,A D *) 

where A r is the r helicity, Ad* is the D* polarization, is the B meson mass, q^ = 
(pb — Pd(*)Y an d Pb £>(*) are the B,D^ meson four-momenta. The three-body phase 
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space <i$3 is written as 



VQ7Q- 

256vr 3 



m 



B 



dq 2 d cos # T 



(14) 



where 



q 2 and m D (,) are the meson masses. Hadronic 



amplitudes in the matrix elements Ai = {D^*Hv(\T-L e ^\B) are defined as 

(D(v D )\cfb\B(v B )) = ^B^E[h+{w){v B + v D y + h-(w)(v B - v D )"}, (15) 
(D*(v D .,e)\crfb\B(y B )) = i y /m B m D .h v (w)£ ttvptr e* l/ (v D *) p (v B ) tT , (16) 
(D*{v* D ,e)\c^ 5 b\B{v B )) = y/m B m D .[h Al {w)(w + l)e*" - h A2 {w){e* ■ v B )v B 

-h A3 (w)(e*.v B )v^}, (17) 

where v B = Pb/tkb, v D (») = p B (*) / m D(*) an d w = v B ■ v B - In the heavy quark limit 
(HQL), the form factors become related to a single universal form factor, the Isgur-Wise 
function £(w) PSIEB]: 

h+(w) = h v (w) = h Al (w) = h As (w) = £(w), 

h_(w) = h A2 {w) = (HQL). (18) 

The form factors, including short-distance and 1/rriQ corrections, are known [T7j. Their 
form factors involve the unknown parameters, which have been analyzed [TBI US]- We 
relate the (pseudo) scalar hadronic amplitudes to the (axial)vector hadronic amplitudes 
by using the equations of motion as 

q„{D\cfb\B) = (m b - m c )(D\cb\B), (19) 
qfl {D*\c^ 5 b\B) = -(m b + m c )(D*\c l5 b\B). (20) 

The other hadronic amplitudes are equal to zero due to parity and time-reversal in- 
variance, i.e., (D\c^b\B) ={D\ cj^'j^b \B) = (D*\cb\B) = 0. See the Appendix for more 
details. 



3 Numerical results 

We evaluate R(D^) and the forward-backward asymmetries as functions of Cs RL on 
B — > D^tPi by using heavy-quark symmetry with short- distance and 1/wiq corrections 



as 



R(D) 
R(D*) 
A FB (D,n) 
A FB (D,p) 
A FB (D,ai) 
A FB (D*,ir) 
A FB (D*,p) 
A FB (D*, ai ) 



0.310(11)] ft, 
0.253(3)] 7T, 
0.54 + 4.0Re(C 
0.32 + 2.4Re(C 
0.25 + 1.9Re(C 
0.28 + 1.3Re(C 
0.092 + 0.79Re(C 
-0.055 + 0.62Re(C~)l / K 



n, 
/n, 
/n, 
/n\ 

I TV 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 
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Figure 2: We fix Im(Cs R h ) = 0. In the left (right) panel, the blue line shows the Cs R . L dependence 
of R(D) (R(D*)), the red line shows the Cs R L dependence of A FB {D, w) (A FB {D* , 7r)), and the light 
blue band corresponds to the measurement of R{D) (R(D*)) at 95% C.L.. 
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Figure 3: The Cs RL dependence of 7r)(red solid lines), A FB {D^*\ p) (green, dashed lines) 

and Afb(D(*> , ai)(gray dotted-dashed lines). 



where 

K = 1 + 9.0Re(C+) + 37|C+| 2 , (29) 
TV = 1 + l.lRe(C-) + 3.9|C~| 2 , (30) 

C + = GeV 2 x ( Cs » + Cs A (31) 

V m b -m c J 

C~ = GeV 2 x ( Cs »- Cs A , (32) 

V m b + m c J 

and m^c are the b, c quark masses. We use the m& and m c in the MS scheme at the mj, 
scale [20] in this paper's figures. A few percent errors due to the measurements and the 
hadronic uncertainties remain. These quantities determine Ke(Cs RL ) and \lm(Cs R ± 
Cs L )\- 

In Figure El we fix lm(Cs RL ) = 0. In the left (right) panel, the blue line shows 
the Cs RL dependence of R(D) (R(D*)), the red line shows the Cs RL dependence of 
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Figure 4: We fix Im(Cs RL ) = 0. In the left panel, the light blue or light red regions show the 99% 
C.L. allowed regions for R(D) or R(D*). In the right panel, four green regions are the 99% C.L. allowed 
regions for R(D) and R(D*). The regions (i), (ii), (iii) and (iv) are classified by Aps as seen in Eq. 
(ED-®. 



Afb(D,tt) (A fb (D* ,tt)), and the light blue band corresponds to the measurement of 
R{D) (R(D*)) at 95% C.L.. In Figure [3J we show the results of the A FB for all modes. 
In Figure HI we fix lm(Cs RL ) = 0. In the left panel, the light blue or light red regions 
show the 99% C.L. allowed regions for R(D) or R(D*). In the right panel, the four green 
regions are the 99% C.L. allowed regions for R(D) and R(D*). This fourfold ambiguity 
cannot be solved by using only R(D^). However, A FB can among discriminate these 
regions. For example, the regions (i), (ii), (iii), and (iv) are classified by A FB (D^*\ ir) 



as 



A FB {D, 7r) > 0, A FB {D\ix) > 0.1 (i), (33) 

Afb{D, 7r) > 0, A FB {D*,n) < 0.1 (ii), (34) 

A FB (D,7t)<0, A FB (D*, 7r) > 0.1 (iii), (35) 

A FB (D,7t)<0, A FB (D*, 7r) < 0.1 (iv). (36) 

Since a tauon decays into a light meson(lepton) with nutrino(s), it is difficult to 
measure the forward-backward asymmetries for tauons and mesons : 

A ^ Di)) = ~ r ( Lwr>) ' (37) 

where 9 T is the angle between the tauon and meson as seen in Figured! However, it 
is not impossible to analyze A FB (D^) by using information about the position where the 
tauon decays. An analysis to reconstruct the tauon would start at the LHCb experiment 
[21] . Then, we evaluate A FB (D^) as functions of Cs RL as 

A FB (D) = [0.358(1)] [1 + 7.1Re(C + )] / TZ, (38) 
A FB (D*) = [-0.065(8)] [1 - 13Re(C - )] / TV. (39) 
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4 Conclusion 



We have studied the decay modes B — > D^rvg with the charged scalar effects, and 
show that it is possible to determine Re(Cs RL ) and \lm(Cs R ±Cs L )\ with the combina- 
tion of the ratios of branching fractions R(D^) and the forward-backward asymmetry 
A fb (D ( -*\tt), A FB (D ( -*\p), and A FB (D ( -*\ai). When considering the effective weak 
Hamiltonian (JBJ), we evaluate R and A FB as functions of Cg RL in Eqs. (!2"Tl) - (f2"8"l) . 

As seen in Figure HI four allowed regions for R(D^) exist. This fourfold ambiguity 
cannot be solved by using only R(D^). However, A FB can discriminate among these 
regions, because the Cs RL dependence of R and A FB are different, as seen in Figure [2J 
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Appendix : Form factors 

In this paper, we use the B — > form factors estimated by the heavy-quark symmetry 
with both short-distance and l/rriq corrections [17]. For B — > Dtu t , we define the 
hadronic amplitudes as 

(D(v D )\cb\B(v B )) = ^b^E(w + l)h s (w), (40) 
(p{v D )\crfb\B{v B )) = V^B^B[h+(w)(v B + v D f + h_(w)(v B - v D )% (41) 

and more, the combinations which appear in the calculations as 

V x {w) = h + (w) - ±ZLh4w), (42) 

S^w) = h + (w) - i±I ^h^w). (43) 
1 — r w + 1 

Vi(w) is parameterized as 

V^w) = Vi(l)[l - Spjz + (51p? - 10)z 2 + (252p 2 - S4)z% (44) 

where z = (y/w + 1 — \/2)/(\/w + 1 + a/2). The parameters Vi(l) and pi have been 
analyzed by the distributions dT(B — > D(e or p)v)/dw, and we use p\ = 1.18±0.06 [T8] . 
The Vi(l) dependence cancel out in the calculations of R{D) and A FB (D). We estimate 
Si(w) as 

Sxiw) = 1.0036[1 - 0.0068(w - 1) 

+ 0.0017(w - l) 2 - 0.0013(w - l^V^w). (45) 

We relate hs(w) to Si(w) by using the equations of motion ffT9l) as 

h s (w) = S^w). 46 

m b - m c 
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For B — > D*tu t , we redefine the hadronic amplitudes as 

(D*(p* D ,e)\c l5 b\B(p B )} = f P (w)(e*-p B ), (47) 

(D*(p* D ,e)\cry»b\B(p B )) = if v {w)e^el{p B + p D *) p (p B - p D *) a , (48) 

(D*(p* D ,e)\cYl 5 b\B(p B )) = f Al (w)f + A»(e* ■ P b)( P b +Pd*Y (49) 

+ f A3 (w)(e* ■p B ){p B -p D *Y- (50) 

We rewrite these form factors to more useful forms as 

f Al (w) = y/m B m D *(w + l)Ai(w), (51) 

- ~2^=k Mw) ' (53) 



Ai(iu), i?i(w), R 2 (w) and Rs(w) are parameterized 



as 



H = Ai(l)[l - Sp 2 ^ + (53p 2 Al - 15)z 2 + (231p 2 Al - 91)z% (55) 

R^w) = i?i(l) - 0.12(w - 1) + 0.05(w - l) 2 , (56) 

R 2 (w) = R 2 (l) + 0.11(™ - 1) - 0.06(u; - l) 2 , (57) 

R 3 ( w ) = R 3 (l) - 0.03(w - 1) + 0.02(w - l) 2 . (58) 

The parameters p\ , -Ri(l) and i^(l) have been analyzed by the distributions dT(B — > 
D*{e or n)v T )/dw, and we use p 2 ^ = 1.214 ±0.035, i?i(l) = 1.401 ± 0.038 and R 2 (l) = 
0.864±0.025 [19]. We estimate R 3 (l) ~ 1.12, and the relation between i2 3 (l) and i2 2 (l) 
as -Rs(l) — -^2(1) + 0.85m£)*/ms from Ref [17]. From the latter relation and experiment 
results, however, ^3(1) are estimated as #3(1) — 1.19 ± 0.03. So, in our calculation we 
estimate as i? 3 (l) = 1.17 ±0.05. Finally, We relate fp(w) to f Al (w), f A2 (w) and f A3 (w) 
by using the equations of motion f[2"0l as 

fpH = ^ [JaM + (m 2 B - m 2 D )f M (w) + q 2 f A3 (w)}. (59) 
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